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Abstract. We present further evidence that gelation is an arrested phase 
separation in attractive colloid-polymer mixtures, based on a method combining 
confocal microscopy experiments with numerical simulations recently established 
in Nature 453, 499 (2008). Our results are independent of the form of the 
interparticle attractive potential, and therefore should apply broadly to any attractive 
particle system with short-ranged, isotropic attractions. We also give additional 
characterization of the gel states in terms of their structure, inhomogeneous character 
and local density. 



1. Introduction 

The addition of non-adsorbing polymers to a colloidal suspension can induce an effective 
depletion attraction between the colloidal particles pQ, leading to a much richer phase 
behavior [21 E] than that normally observed for atomic or molecular systems [I]. Colloid- 
polymer mixtures are therefore often used as model systems of complex fluids, which 
are normally encountered in technological applications such as food and personal care 
products, where the attractive interactions may significantly impact product stability 
Micron-sized constituent colloids are large enough to be imaged with light, while still 
small enough to have their dynamics driven by fc^T, the thermal energy, where T is the 
temperature and fee, the Boltzmann constant. Both the range and magnitude of the 
attractive interactions can by finely tuned by adjusting the radius r p and concentration c p 
of the added polymers. With sufficiently high c p , these mixtures form gels, i.e. networks 
of particles that span across the volume in a disordered arrested structure, which can 
sustain shear stresses, despite the small packing fraction occupied by the colloidal 
particles [5j [6j [7] . In the limit of strong attractions and very small 0, gelation is driven by 
the diffusion-limited cluster aggregation (DLCA) of fractal clusters of colloidal particles, 
a mechanism well-studied in the past [HI El HO]- In the limit where <fi is large, particles 
arrest to form colloidal glasses[llj, which can be of two fundamental types: repulsive, 
driven only by packing and geometrical frustration [T2|, and attractive, when particle 
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bonds become dominant, at sufficiently high attraction strength[13j Q3]. A simple 
extension of attractive glasses to lower and lower densities is not necessarily a good 
explanation for gel formation, due to the intervening of liquid-gas (colloid rich-colloid 
poor) phase separation. Several different mechanisms have been proposed for the gel 
transition at low- and intermediate-^ and high interaction strengths |7j. These include 
percolation [151 [16], mode-coupling theory (MCT) attractive glass [T7]. cluster-MCT[T8] . 
fluid-crystal transition[19j, spinodal-induced gelationJ^Ul EH [22] , viscoelastic phase 
separation [23]. Stimulated by a series of numerical studiesj2H [25j [26] and by analogies 
with studies of atomic fluids [27] and polymer solutions [2H], the hypothesis of arrested 
spinodal decomposition has recently received a particular attention in experimental 
studies [29], [301 EI]- One of the greatest challenges in reaching an unambiguous 
interpretation is to find a way to precisely quantify the attractive interaction potential, 
which would allow a quantitative connection to be forged between observed experimental 
phenomena and a well-characterized thermodynamic phase or state. 

Thus far, the mapping of the experimental c p into a system-independent 
thermodynamic parameter has been based on assumptions of a specific shape for 
the interparticle potential V(r), such as the Asakura-Oosawa (AO) model[l] typically 
employed to describe colloid-polymer mixtures. Within this model, colloids (with 
radius a) behave as hard spheres, while polymers forming coils in solution have a dual 
nature: with respect other polymers, they are an ideal gas in which polymer coils can 
interpenetrate; with respect to the colloids, polymers behave as hard spheres of radius 
r p , which cannot interpenetrate the colloidal spheres. With these simple assumptions, 
an effective colloid-colloid attractive potential is derived, 



where II is the osmotic pressure, V the volume and F(x) a well defined functionp] which 
vanishes for x > 1. The maximum strength of attraction occurs when two colloidal 
spheres touch, and this attraction is denoted by the parameter U = Vao( t — o)/kBT. 
To facilitate comparison between different experimental conditions, the concentration- 
dependent radius of the polymer is typically divided by the colloid radius to yield 
a dimensionless size ratio £ = r p /a. A universal gelation state diagram may then 
be parameterized by 0, £ and U, yielding a way to compare different experiments 
and theories [52]. However, while <fi and £ are unambiguously defined in experiment, 
simulation and theory, U/ksT depends entirely on the shape of the interaction potential. 
This model-dependence of U leaves its quantification vulnerable to the uncertainties in 
the applicability of a given theory. 

Recently Noro and Frenkel[33] have established a generalized law of corresponding 
states for systems interacting with short-range attractive potentials (£ < 0.1), in terms 
only of the normalized second virial coefficient B\ = (3/8a 3 ) / °°(1 — e~ v ^ r ^ kBT )r 2 dr. 
For a given value of B%, any system interacting with any given short-ranged attractive 
potential (e.g. AO, square well or generalized Lennard- Jones) shares the same 
thermodynamic properties and structure. Hence, either a direct measurement of B\ 




(1) 
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or a connection of the experimental attraction strength (e.g. expressing c p in terms of 
B2) would allow a universal description of the system. For example, one recent work on 
lysozyme has used the Noro-Frenkel law to express the crystallization transition under 
various solvent conditions in terms of a universal -B^IM]- Another possibility often 
exploited in the past is to rely on the well-known Baxter model[35], which corresponds 
to the sticky hard sphere limit: the attraction range becomes infinitesimal while the 
attraction strength becomes infinite in such a way that the second virial coefficient 
B2 remains finite. This model is theoretically convenient because it is amenable to 
analytic treatment within Percus-Yevick (PY) closure [36J and, consequently, it has been 
used for a long time in comparison to experiments |37] . However, the thermodynamic 
inconsistency of the PY solution has often led to an inaccurate location of the gel line 
with respect to the phase-separation line in the thermodynamic phase diagram [38J. 
Only a recent accurate location of the gas-liquid phase separation, which was obtained 
numerically by Miller and Frenkel|39j, has opened the way for a more careful use of the 
Baxter model, via the B\ scaling, in relation to experiments where samples were shown 
to reach equilibrium over the course of several months [30J. 

Only a truly model-independent measure of the attractive potential, that could 
be applied to a variety of attractive colloidal systems, would facilitate quantitative 
comparison between theory and experiments, and help to elucidate the mechanisms 
driving the onset of gelation. To address this issue, in a recent paper [3D], we have 
introduced a novel method for generic short-range isotropic attractive interactions - 
independent of any specific potential shape or model — that is capable of quantitatively 
mapping the experimental observations of gel formation in colloid-polymer mixtures 
onto the corresponding thermodynamic phase diagram. Through the comparison of the 
cluster mass distributions, measured directly at the single-particle level using confocal 
microscopy and in simulations of various short-ranged attractive potentials, we have 
established a one-to-one relationship between the experimental c p and B%, building on 
the Noro and Frenkel law|33j. In doing so, we provide unambiguous evidence that 
gelation takes place precisely where spinodal decomposition occurs, evidencing a strong 
link between thermodynamic behaviour and dynamic arrest. The gel is triggered from 
the initial separation into two phases, of which the densest may become locally arrested 
due to an intervening attractive glass transition [7], and hence kinetics only intervenes 
after thermodynamics has provided the necessary density fluctuations. 

In the present manuscript, we report additional analysis of the experimental system 
in both the fluid and the gel phase to complement the results previously reported in[40j. 
These results provide further evidence that gelation in short-ranged attractive colloid- 
polymer mixtures is an arrested phase separation. 

2. Description of the experimental system 

Sterically-stabilized colloidal spheres of polymethylmethacrylate (PMMA) with DilCig 
fluorescent dye are suspended in a 5:1 (by mass) solvent mixture of bromocyclohexane 
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(CXB, Aldrich) and decahydronaphthalene (DHN, Aldrich) for several months, following 
procedures already described in [32j HOj HI]. The PMMA particles have radius 
a = 560 ± 10 nm, determined with dynamic light scattering. Tetrabutylammonium 
chloride (TBAC, Fluke) is added until saturation (~4 mM) to screen any residual long- 
range Coulombic repulsion. The colloidal suspension is carefully buoyancy-matched [40J. 
Two sets of depletant polystyrene (PS) are used, with molecular weights respectively of 
Mvi/=69.2 kDa and M^=681 kDa. From Zimm plots of static light scattering data, the 
radii of gyration of the two polymers are determined to be 10.0 and 33.0 nm respectively, 
yielding £ = 0.018 and £ = 0.059 in the two cases, always well within the short-range 
attraction limit where the Noro-Frenkel law is valid. The overlap concentrations, defined 
as c* = 3Mw / i^TpNA, are 27.2 and 7.5 mg/mL, respectively for the two sets of polymers, 
with Na being the Avogadro's number. In all cases, the raw polymer concentrations 
are directly measured as the mass ratio of mg PS per g of total sample mass, then 
converted into a ^-dependent free- volume c p (mg/ml)[3]. Samples at varying polymer 
concentration c p are generated while preserving constant 0, although a small decrease is 
observed with increasing polymer concentration [JO]; this may be due to larger clusters 
sticking to the walls of the macroscopic sample chamber. 

Following our previously-reported imaging protocol [321 HUJ, HI], 3D stacks of 181 
images are collected, each an 8-bit TIFF of 1000 x 1000 pixels, at 10 frames per second. 
Each image stack covers a volume of 60 x 60 x 60 /im 3 , taken from the center of the 
sample at least 20/zm away from any sample chamber surface to minimize edge effects. 
Each fluid sample is homogenized and equilibrated for four hours, then 26 independent 
3D image stacks are collected with an automated confocal microscope[41J. Each gel 
sample is homogenized and then observed immediately, with 3D stacks of the same 
sample volume collected every 50 seconds for the first 5000 seconds, then every 1000 
seconds for the subsequent 100,000 seconds. 

Using confocal microscopy, we observe and individually locate each colloidal particle 
in the sample imaging volume in 3D; a typical gel is shown in Fig. [TJ Most of the colloidal 
particles are connected in a spanning network, while at all times a residual fraction of 
particles is present in monomers and smaller clusters that freely diffuse throughout the 
voids of the network. We also observe a continuous exchange between these unconnected 
particles and the gel network, due to the finite attraction strength of order a few ksT . 
However, once a percolating gel network forms (which happens generally within a few 
thousand seconds), it persists for all experimentally observed time — it is arrested. 

Locating each particle individually with confocal microscopy also allows more 
detailed structural analysis, such as a direct calculation of the static structure factor, 
defined as S(q) = (| J2iLi e iq ' rj | 2 )/iV. For fluid samples, we average over the 26 
independent configurations. For gel samples, we follow a single configuration over time. 
We calculate S(q) for a cubic box of size 50/im, centered within the sample imaging 
volume, in order to minimize edge effects, which if present, will only affect the range 
2qa < 0.2. 

This precise way to determine the static structure factor allows a very stringent 
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Figure 1. 3d-Reconstruction of the gel for (f> = 0.14, £ = 0.059 and c p = 1.80 mg/ml. 
Isolated clusters are shown in different colors for different cluster sizes, while monomers 
are depicted as transparent (and smaller) spheres. 




Figure 2. For <f> — 0.16 and c p = 0, (a) static structure factor and (b) direct 
correlation function, calculated using the confocal microscopy coordinates (circles), 
compared to the PY calculations for hard spheres at the same <f> (full lines). 



test on how close our colloids are to being a model hard-sphere system. In FigMa), 
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we plot S(q) for an experimental sample with = 0.16 and c p = 0, together with PY 
theoretical predictions at the same <fi. The comparison also permits an estimate of the 
particle radius which best reproduces the superposition of the peak positions, which in 
this case turns out to be ~ 600nm, quite close to the one measured with dynamic light 
scattering. We observe a remarkable agreement of S(q) with the theoretical predictions 
essentially over all length-scales, which demonstrates that, in the absence of polymer, 
our colloids behave as hard spheres. Similar results are found for the other investigated 
values. We also examine the direct correlation function g(r), the Fourier transform 
of ^(^[SS]- Experimental data and the PY prediction are shown for the same sample 
in Fig [2^b). Because these fluid particles diffuse while being imaged, their motion 
provides a fundamental physical limit on how precisely their centers can be located. Our 
uncertainty from the image-processing software is well under 50 nm in all directions [42J. 
This positional uncertainty smears out the shape of g(r), especially close to the contact 
peak, resulting in an apparent worse agreement with the PY solution. We therefore 
conclude that the S(q), whose signal is maximum around qa ~ n, is a better observable 
for direct quantitative comparison between confocal images and theory/simulations. 

3. Locating the gel transition onto the thermodynamic phase diagram 

To locate the gel transition on the thermodynamic phase diagram, we have employed 
a method based on the comparison of experimentally-determined cluster mass 
distributions with those obtained from simulations in the fluid phase up to gel point [40J. 
For each experimental sample (with a well-defined c p ), its cluster mass distribution 
was determined, after defining as connected any pair of particle whose distance is less 
than 1.3/im. The choice of this value, which is larger than the range of attraction, 
allows us to minimize spurious effects introduced by the aforementioned particle position 
uncertainties. Simulations were run at the same value of <fi and several values of B^; in 
order to identify the optimal B\ value, the cluster mass distribution that best matches 
the experimental data is identified via a least-squares minimization. Connected particles 
were defined with the same procedure in both experiments and simulations. There are 
no free parameters involved in this matching procedure, which uniquely associates each 
experimental c p with a theoretical value of B\. Because the experimental conditions are 
in the short-ranged attractive limit, we can in principle use any short-ranged attractive 
potential in the simulations. By using several potential shapes (AO, square well and 
generalized Lennard- Jones) — all of which yield quantitatively identical results [4TJj — we 
provide stringent reconfirmation of the validity of the Noro-Frenkel law. 

The results of this c v -B\ mapping procedure are presented in Fig. [3] for <fi = 
0.045 and the two studied sets of polymers described above. Horizontal lines mark 
the experimentally-observed concentration threshold c 9 p el denoting the transition from 
fluid to gel for each polymer. The vertical line denotes the boundary of spinodal 
decomposition expressed as a value of B%, which depends only on </>. This boundary is 
estimated both within simulations, through the monitoring of the time-dependence of 
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Figure 3. Dependence of polymer concentration on the thermodynamic universal 
parameter B% for <fi — 0.045 and two sets of polymers, with £ = 0.018 and £ = 
0.059. The experimental gel transition is observed to happen exactly at the onset of 
thermodynamic phase separation. The dashed line is the AO prediction for £ = 0.059. 



the energy and of the small angle static structure factor [21], and via the PY energy 
route [SB], which has been shown to be quite accurate for predicting the spinodal 
temperature [39]; both values coincide within the error bars. Without exception, the data 
reported in Fig. [3j as well as those for other studied <p (not shown here) [ID] , demonstrate 
that the gel transition takes place right at the onset of spinodal decomposition. (Note 
that the B\ scale is rather expanded when compared to a temperature scale, so that 
the observed distance from the last fluid sample to the gel boundary is very small, 
and comparable to the experimental error in c p .) These data suggest a mechanism for 
gelation of short-ranged attractive particles: the system first enters the unstable region, 
where density fluctuations drive the formation of locally-dense liquid droplets; then 
these denser regions arrest probably through an intervening attractive glass transition. 

Fig. [3] provides the 'true' relationship between polymer concentration and 
thermodynamic parameters: given the correct potential model to describe the system, 
its B\ should superimpose onto these points. However, since B\ is an integrated quantity 
of the potential in terms of particle distance, the original potential cannot be recovered 
easily from B\ alone. However, if the simple AO formula (Eq{T]) is used to determine 
B\ as a function of c v (dashed lines in Fig. [3]), the resulting relation is rather different 
from what the cluster mass distribution mapping predicts. To reach any given B%, 
the AO model predicts a much smaller value of c p than that observed in experiment, 
demonstrating that it does not quantitatively describe the potential for our experimental 
systems. The discrepancy, as expected, becomes larger and larger for increasing polymer 
concentrations. 
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4. Characterization of the Gel States 

We now examine some features of the observed gel states. For each studied 0, we have 
monitored the evolution over time of several gel states at different c p . In Ref.[40j we have 
reported S(q), obtained from the particle coordinates, across the fluid-gel transition. We 
observed that while the fluid, even at the highest c p , retains a low value of S(q — > 0), a 
big jump of about two orders of magnitude in the low-g peak of S(q) is found as soon 
as the system crosses to a gel state. Moreover, we followed the evolution in time of 
this spinodal peak, finding an initial slow coarsening process (compatible with a power 
law t -1 / 6 ), later interrupted by arrest, i.e. at long times the structure of the system 
reaches a final state and remains subsequently unchanged [40J. These results are in good 
agreement with numerical simulation studies in the arrested spinodal region[24, 25J. 

We here complement these results with the additional evidence that the arrested 
structure crucially depends on the location of the studied state point within the two- 
phase region. Indeed, we plot in Fig|4](a) and (b) the evolution of S(q) with time for 
= 0.16 and two different c p values. On the left panel, the data correspond to a system 
which is just inside the gel boundary, i.e. a (shallow) quench just within the spinodal 
region. In this case, the low-g peak of S(q) undergoes quite a significant coarsening 
process, growing in amplitude by almost two orders of magnitude, before eventually 
arrest takes place and the growth stops at long times. On the right panel, we report 
data corresponding to a much deeper quench inside the spinodal region. Here, the 
coarsening process almost does not proceed at all, and the structure stops changing 
quite rapidly, without even developing a large compressibility. Also the peak position 
does not move significantly with time. These results show that the resulting gels formed 
by arrested spinodal can differ very much in their final structure, and in particular in 
their mesh size and compressibility, depending on the chosen quench path inside the 
two-phase region. 

We also monitor the evolution of S(q) using the scaling form for spinodal-driven 
aggregation derived by Furukawa[4"3"]. 

S{q/qi,t)~q x {t)- d F{q/ qi ), (2) 

which was also generalized to the case of fractal clusters undergoing DLCA 
by Carpineti and GigliofTU]. Here as typical length scale we use qi(t) = 
(f Q 9c S(q,t)qdqj / (J 9c S(q,t)dq) where q c = 6a is an appropriate cut-off value, which 
provides the first moment of the wave- vector in the low-g peak of S(q) to reduce 
numerical noise, while F(q/qi) is a time-independent scaling function and d is the 
spatial dimension. For fractal aggregation this is found to be replaced by the fractal 
dimension of the clusters[TU]. Within our time window of observation of the gel states, 
three different time intervals can be identified: an early stage (which we observe only 
partially, as our first observation time is equal to 100s), an intermediate regime where 
the Furukawa scaling should hold; and, finally, the late stage or arrested regime, where 
the system structure does not evolve any longer. In order to visualize the intermediate 
regime and to demonstrate the scaling of Eqj2]with d = 3, we focus on state points not 
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Figure 4. (a) Shallow versus (b) deep quench evolution in time of S(q) for <f> = 0.16. 
(c) Scaling of S(q) as in [10 for the shallow quench case reported in (a) and t > 3000s. 
The inset shows the time dependence of q\ for the same state point; at short times the 
reduction of the characteristic length scale is compatible with t -1 / 6 , later followed by 
a slower decrease. Symbols are consistent throughout all panels. 

too far from the spinodal boundary, where we observe a significant evolution of S(q) 
over time. We find that our data, under these shallow quench conditions, do scale onto 
a single master curve at intermediate times, well before arrest takes place, as predicted 
by the scaling. Trivially, data from later times also scale onto the same master curve. 
As an example we show in Figj4](c) the scaled data corresponding to the shallow quench 
of Fig|4](a). For this particular state point, a significant evolution of S(q) over time 
is observed, as it is evident from the time-dependence of q\ reported in the inset. For 
t > 3000s, which interestingly appears to be very close to the time where q\ departs 
from the t~ 1 ^ law towards a slower form (until saturation, not observed just for this 
particular sample), the scaling holds at all subsequent times. For all other studied state 
points, arrest takes place much earlier so that the scaling can be demonstrated only in 
a narrower regime. In all cases, the scaling does not providence evidence of a fractal 
structure of the clusters. 

Another measure of the inhomogeneous character of the arrested states is the 
so-called demixing (or inhomogeneous) parameter \EvlII], which can be evaluated by 
dividing the system into n 3 boxes and measuring the local density in each box. Hence, 

n 3 

* n = ^|p fc -p| 2 , (3) 

k=0 

where pk is the density in the fc-th box and p = 3(f)/ (4na 3 ) is the average particle 
number density. In a fluid or homogeneous state, this parameter averages to zero, 
while as soon as inhomogeneities appear, it grows. A comparison of \&„ for different 
attraction strengths can provide information on how the degree of inhomogeneity varies 
with different quench paths. We show results for \& 4 — calculated by dividing the imaged 
volume of (52pm) 3 into 4 3 = 64 boxes of 13/xm edge length — versus time, in Fig|5[a) for 
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<p = 0.045 and all studied quench depths. We also repeated our analysis for a smaller box 
edge (of 4/xm), finding the same qualitative results. The evolution in time of ^4 for all gel 
samples displays an initial growth followed by a plateau region, similar to what observed 
for the maximum height of the low-g peak of S(q) , named S MAX (q), for the same state 
points, reported in Figj^b). In the inset of Fig|5^a) the average (for fluid states) and 
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Figure 5. (a): Demixing parameter ^4 evolution with time for the gel states with 
4> = 0.045. Inset: the fluid to gel transition in terms of ^4; (b)S MAX (q), the maximum 
height of S(q), evolution with time for the same state points. Identical lines and 
symbols are used in both panels. 

long-time limiting value (for gel states) of \&4 is reported, in order to show, once again, 
the crossing of the spinodal region in correspondence of the gel transition. It is intructive 
to notice that, apart from the gel state very close to the spinodal boundary (i.e. with 
c p = 3.31mg/ml), whose spinodal decomposition process seems to proceed further than 
the others before finally arresting, all other gel states, corresponding to deeper quenches, 
seem to converge toward a similar value of the inhomogeneity parameter. Similar results 
are also found for <fi = 0.16. 

Finally, we examine the local density (or <fi) distribution for several values of the box 
size. This approach is sometimes implemented in simulations to establish the presence 
of a two-phase region[l5]: if two distinct peaks appear, then two phases of different 
densities can be properly identified. The observation of the two peaks is strongly 
connected to the choice of a box size comparable to the size of the spatial density 
fluctuations, and requires a significant evolution of the phase separation process to be 
observed[13]. Fig. [6] shows the evolution in time of P((f>) for two different box sizes, 
respectively L = 3//m and L = 5//m for = 0.045 (corresponding to an average number 
of particle per box of respectively 1.5 and 7.1). Results show that the system starts from 
a rather homogeneous condition (since P(4>) on the L = 5/im scale shows a maximum 
around the average 0), then evolves in the direction of increasing the density fluctuations 
(since states with both small and large packing become preferentially probed) until arrest 
is reached. On the small L = 3 scale, the local density increases even more than a factor 
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15, reaching values comparable to a dense packing of spheres (0 > 0.6), consistent 
with the estimates of the internal gel density derived in Ref. |40| . Smaller <p values are 
probed on the L = 5 scale, suggesting that the characteristic mesh of the aggregate is 
comparable to the smaller box size. We also note that, in the time window preceding 
dynamic arrest, we do not observe the formation of a bimodal distribution, suggesting 
that gel formation may precede true bulk phase separation. 




Figure 6. Evolution in time of local packing fraction probability P(4>), obtained 
measuring the local <f> in small boxes of sizes 3 (left) and 5 (right) /im), for the state 
point <j) — 0.045 and c p = 3.31 mg/ml. 



5. Conclusions 

In this manuscript, we have provided additional evidence of the mechanism at the 
origin of the gel transition in short-ranged attractive colloid-polymer mixtures. We 
have firstly shown that the experimental colloidal system is well characterized: without 
added polymers, it behaves as a system of nearly-perfect hard spheres. Then, comparing 
cluster distributions in experiments and simulations has allowed us to precisely map 
the experimental gel transition on the thermodynamic phase diagram, quantifying 
the attraction strength in terms of the normalized second virial coefficient B\. This 
universal parameter contains all the necessary information for short-ranged isotropic 
attractive potentials to entirely determine the structure and connectivity properties of 
the system. In this way, our procedure is independent of any specific choice of interaction 
potential model and is relevant to any short-ranged isotropic attractive system. Our 
findings unambiguously show that thermodynamic phase separation triggers the kinetic 
arrest [ID] . 

A growing consensus has recently favored the arrested phase separation mechanism 
for gelation in these short-ranged attractive systems. A number of simulations 
studies [2H [251 126] and various experimental colloid investigations — including colloids 
with depletion interactions [23 EDJ as well as mixtures of oppositely charged 



Gelation as arrested phase separation in short-ranged attractive colloid-polymer mixturesYl 

colloidsjl6] — all support this picture, which may also be relevant for atomic glass- 
formers^]. Moreover, our results confirm important conclusions drawn from data on 
attractive lysozyme solutions, recently reported by Cardinaux et al. [31]. Both in our 
colloidal system and in the protein solution case, gelation is caused by the arrest of the 
dense phase induced by the liquid-gas spinodal decomposition. Yet differences between 
the two systems appear to arise in the estimated value of the gel packing fraction; 
this packing fraction appears to be quench-depth independent as evaluated in Ref. [10] 
and further confirmed by the present analysis, while dramatically depending on quench 
depth in lysozyme, where the density of the dense phase shows a marked decrease 
with increasing attraction strength [3 1 j . However, several factors may contribute to 
these reported differences. First, differences in the interaction potential are present 
between lysozyme and colloid-polymer mixtures. In lysozyme, electrostatic repulsion[47j 
can play an important role and, in principle, anisotropy effects (due to shape or 
patchyness) could be relevant. Indeed, one observes a rather low-0 position of the 
critical point (0 C ~ 0.18[3l]) as compared to standard short-ranged attractive potentials 
(0 C ~ 0.27(18]). Second, the methods used to establish the local gel volume fraction are 
rather different in the two cases. While we have used a microscopic method, based 
on single-particle location, to estimate the internal volume fraction of the gel |40j, 
Cardinaux et al. determine the gel volume fraction by centrifuging the lysozyme solution 
after the phase separation was arrested, decanting the supernatant, and measuring the 
resulting macroscopic volume fraction. Whether these two metrics would yield the same 
results, if they could both be performed on the same sample, it is presently unknown. A 
further complication arises by the fact that single 3.4 nm lysozyme proteins cannot be 
individually resolved. Finally, the finite size of the imaged volume in our experiment, as 
compared to the macroscopic volume (compared to particle size) accessed in [31] , could 
also be used as an explanation for the discrepancy, although we performed additional 
checks for larger samples of millimetric size confirming our previous findings. 

Another interesting question relates to why we do not observe any state 
point undergoing a complete gas-liquid phase separation within our study. In all 
experimentally studied cases for samples that begin to phase separate, we find that 
a gel forms at long times. Indeed, for such short-ranged attractive systems at the values 
of <p that we study, the difference in attraction strength that would allow to explore 
the liquid region is very small and roughly comparable to the experimental error in c p . 
However, we do observe some state points, very close to the spinodal boundary where 
arrest only arises at very long times (see Fig|4]). Also, it appears that our studied state 
points always refer to quenches inside the spinodal, where the liquid droplets rapidly 
coalesce and form an arrested spanning network. Another interesting case would be 
to look at quenches inside the region at low delimited between binodal and spinodal 
lines, where unconnected clusters may become arrested, forming the so-called glassy 
'beads' [HI i9]. 

Finally, we discussed in some detail the properties of the observed gel states. We 
have provided evidence, from complementary types of analysis, that the resulting states 
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are rather inhomogeneous, that their final structure depends on the quench path, and 
that locally the internal gel volume fraction is very high |40j. Kinetic arrest intervenes 
before phase separation has proceeded enough to allow the visualization of two distinct 
peaks in the local density distribution [45], and the system becomes trapped in a highly 
heterogeneous state. Future work will be devoted to attempts to validate our method 
at larger packing fractions and to test its limit of validity with increasing range of 
attraction. 
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